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The antiferromagnetic Heisenberg two-leg ladder in the presence of frustration and an external 
magnetic field is a system that is characterized by two sorts of quantum criticalities, not only one. 
One criticality is the consequence of intrinsic frustration, and the other one is a result of the external 
magnetic field. So the behaviour of each of them in the presence of the other deserves to be studied. 
Using the Jordan- Wigner transformation in dimensions higher than one and bond-mean-field theory 
we examine the interplay between the field-induced and frustration-induced quantum criticalities 
in this system. The present work could constitute a prototype for those systems showing multiple, 
perhaps sometimes competing, quantum criticalities. We calculate several physical quantities like 
the magnetization and spin susceptibility as functions of field and temperature. 
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FIG. 1: The two- leg ladder showing the couplings along the chains, rungs, and diagonals is 
displayed. 



I. INTRODUCTION 



In low-dimensional antiferromagnetic (AF) Heisenberg spin ladders, the mixing of mag- 
netic frustration, external magnetic fields, and thermal fluctuations together with quantum 
fluctuations is a recipe for interesting non-conventional and exotic behaviours. The two- 
leg frustrated ladder is an example, which is characterized by three distinct non-magnetic 
quantum spin liquid states; the Neel-type (N-type) state, ferromagnetic-type rung (R- 
type) state, and ferromagnetic-type chain (F-type) state. These states are characterized 
by ferromagnetic spin arrangements along the diagonals, rungs, or chains respectively. In 
a recent work, we studied the effect of thermal fluctuations on these states and found that 
the quantum phase transitions between the R-type and N-type states on one hand and 
between the R-type and F-type states on the other hand evolve with temperature into 
classical phase transitions between these disordered states^. Numerical studies^ have 
shown that in the presence of a magnetic field, frustration is able to stabilize a magneti- 
zation plateau at half the saturation magnetization. In terms of the Jordan- Wigner (JW) 
fermions 4 , this would mean that a plateau in the magnetization appears because of the 
appearance of an energy gap between two (excitation) energy bands. In this paper we 
will check this claim within the analytical bond-mean-field theory (BMFT^^i&i. We 
will discuss the interplay between the quantum criticality induced by frustration and that 
produced by magnetic field. Indeed, magnetic field has been shown recently to induce 
quantum criticality between two distinct states, one characterized by spin bond "order" 
and the other one by no spin bond order, in the Heisenberg chain and non-frustrated two- 
leg ladder—. Note that spin bond order does not imply any kind of magnetic long range 
order; i.e., all the states are magnetically disordered in agreement with the Mermi- Wagner 
theorem^. 

In section |TT] the BMFT is briefly reviewed. In section Mil the magnetization, phase 
diagram, mean-field parameters, energy spectra, and susceptibility are examined. In sec- 
tion [IV] conclusions are drawn. 



II. METHOD 



We apply BMFT following the method described in Ref. 7 . The Hamiltonian for the 
spin-i two- leg ladder with diagonal interactions in a magnetic field is written as 

N 2 N N N 2 

H = J Sij-Si + ij+J± S^i-Si_2+Jx / XS^i-Sj-t-i^+Sj-t-i^i-Sj^)— h Sfj, 

i j—1 i i i j—1 

(1) 

where J is the coupling along the chains, J± the transverse coupling, and J x the coupling 
along the diagonals as seen in Fig. Q] The index i labels the position of the spins along 
the two chains, each of which has N sites. The first term represents the interactions of 
nearest-neighboring spins along the chains (legs) of the ladder, the second term represents 
the interactions of the spins along the rungs, and the third term sums the interactions 
along the diagonals. As usual, S,,j is the spin operator. Here h = g^BB with B being the 
magnetic field, g the Lande factor, and fj,B the Bohr magneton. 
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The JW transformation^ for the two-leg Heisenberg ladder is defined as£ 

= °i,j e ^*' 3 > = n i,j ~ 1/2) n i,j — c iJ c i,ji 

i-1 2 i-1 2 

^^^EE™^' ^> 2 = ^EE^-f ( 2 ) 

d=0/=l d=0/=l 

Here i and j are the coordinates along the chain and rung directions, respectively. The 
phases <pi.j are chosen so that all the spin commutation relations are preserved. The c\j 
operator creates a spinless fermion at site while Cjj annihilates one, and n^j is the 

occupation number operator. Using p]). the Hamiltonian ([1]) becomes 

N N 

H =2 E(4i el ™ l ' 2c *+M + ct 2 e^+"c i+1)2 + H.c.) + -i ^(ct lCi , 2 + H.c.) 

i i 
N N 



^^jy^+^W* + c t +iiC . 2 +H.c.) + - \)^+^ \) 

i i j=l 

N 1 1 N 1 1 

- ^)Oi,2 - + ElK.l - ^)Oi+l,2 ™ ~) 



1 1 ^ 2 1 

After applying the JW transformation the Ising terms are decoupled using the Hartree- 
Fock approximation, which neglects fluctuations around the mean field points; (O — 
(0))(0' - (C)) ~ 0, where O and O' are quadratic in d and c 5 . To apply BMFT 
we introduce three mean-field bond parameters; Q in the longitudinal direction, P in the 
transverse direction, and P' along the diagonal. These can be interpreted as effective hop- 
ping energies for the JW fermions in the longitudinal, transverse and diagonal directions, 
respectively^: 

Q = (ci,jcl +1J ), P={cij4j+i), P'=(c t+ i.jcl + i)- (4) 

We choose to place an alternating phase of ir along the chains so that the phase per 
plaqucttc is 7r 10 . This configuration is used to get rid of the phase terms in the Hamiltonian. 
We also set Qij = Qe l ^ i ' i where Q is site independent^. Here <I>j ; j is the phase of the 
bond along the chain such that <f> = ir or 0. This is necessary in order to recover the proper 
result in the limit J x and Jj_ becoming zero, in which we get a result comparable to that 
of des Cloiseaux and Pearson^ for the spin excitation spectrum for a single Heisenberg 
chain, E(k) = | J |sin k\ . 

The magnetization per site is M z = (Sfj) = (riij) — \ . Once a magnetic field is applied 
the magnetization must be included in the decoupling of the Ising terms. The latter must 
be decoupled as many ways as is physically acceptible in order to account for all effects. 
So in addition to the previous two waysi it is now decoupled in the following third way: 

•AS; ,.S';. , , « JM z cj +lj c l+1 ^ + Jc\ ]C ^M z - JM 2 Z - JM Z . (5) 

The mean-field Hamiltonian becomes 

H = J2 *! W + 2N Q 2 + NJ ± P 2 + Nh- N{2J + 2J X + J ± )M Z (M Z + 1). (6) 

k 

As mentioned in the introduction, when J x <C Jj_, the system adopts the N-type state 
with ferromagnetic spin arrangements along the diagonals. When J x 3> J± the system 
adopts the R-type state. In this case, the AF spin arrangements shift to the diagonals, 
and spins on the rungs are forced to adopt a ferromagnetic arrangement. When J x ^> J 
and Jj_ 3> J, the system adopts the F-type state, where the spins on the chains adopt a 
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FIG. 2: In the left panel, the three possible ground states of the system in the Ising limit are 
drawn. In the right panel, the labeling of sublattices corresponding to the short-range spin orders 
that replace the long-range ones in the Heisenberg limit are shown. 



ferromagnetic arrangement and an AF one on the rungs. We use the spin arrangements on 
the right panel in Fig. [2J where the lattice is subdivided into two sublattices, as a starting 
point. The Hamiltonian density for the N-type state is given by 



-hi 



iJisinfc J x iCosfc 



Jj 



H = 



-hi 



1 \ 



ty- J x i cos k 



iJisink 

—hi i<J\smk 
Jxi cosfc — iJisinfc —h' ) 



J x i cos k ^ x 

V i±x 

V 2 



(7) 



where 



Ji = J(1 + 2Q), 

J±i = J±(1 + 2P), 

Jxl = Jx(l + 2P'), 

hi = h-(2J + 2J x +J±_)M Z 

Diagonalizing TL yields the following four eigenenergies 
E Nl {k) 
E N2 {k) 
E N3 {k) 
E N4 {k) 



(8) 



h'+ J x i cosfc + \j Jl sin 2 k + 



J2 

hi— Jxi cosfc + \j J 2 sin 2 k + 



-h'+ Jxi cosfc 
-hi— Jxi cosfc 



J2 sin 2 k+ Jll 

1 4 



(9) 



Similarly, the eigenenergies for the F-type state and for the R-type state are given respec- 
tively by 



J 2 

E F {k) = -ti± Jicosfci t/j 2 lS in 2 fc + - J - L 



J±i 



E B {k) = -h'±^-± ^/J 2 sin 2 k + J 2 ! cos 2 k. 



(10) 



The free energy per site is 



F = JQ 2 + ^ + JxP' 2 + ^-(J + J,+^)M z (M z + l)-^J2Y,Hl + e^% 

k p=l 

(11) 
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FIG. 3: (a) The free energies are plotted as functions of magnetic field for J = 0.5 J±, and 
J x = 0AJ±. (b) The field-induced transitions are shown by zooming in on fields near the critical 
fields. 



where k = N, R, or F depending on the state considered. E Kp (k) designates one of the 
four eigenenergies in state n. The magnetization can be calculated using M z — —^jt, 
which gives 



M, 



(12) 



k p=l 



This adds another equation to the set of self-consistent equation satisfied by the spin bond 
parameters Q, P, and J", which arc obtained by minimizing the free energy: 



Q = - 



i 4 
N.T EE 



p = - 



dE Kp {k) 
8NJ ^ OQ 

k p=l 

1 \ - \ - dE Kp (k) 



4NJ ± 



k p=l 
4 



9P 



n_F[£; Kp (fc)], 
n F [E Kp (k)}, 



P' = — 



JrEE^i^-'M. 



8iVJ 



(13) 



The equations in (|13p are solved numerically in order to obtain the free energies, magne- 
tization, bond parameters, and uniform spin susceptibility as functions of magnetic field 
and temperature. 



III. RESULTS 



A. Free energies 



The free energies are compared as seen in Fig. [3] to determine the state with the lowest 
free energy, which is in turn used to obtain the parameters and magnetization. As Fig[3jD 
shows for the given coupling constants set, a first transition occurs from the N-type state 
to the F-type state, then another one from the F-type state to the R-type state. Note that 
this couple of transitions was not mentioned in Ref. 2 , which used the Lanczos method. 
Our approach works well for analyzing the effect of magnetic field when J x = C£, and 
when h = and J x ^ 1 . Including both J x and h may be pushing the the current 
mean- field approach to its limits. In the strong field regime all free energies are equal 
because the magnetization has saturated and all AF fluctuations have disappeared; i.e., 
even short-range AF correlations are absent in the saturated ferromagnetic state. 
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FIG. 4: (a) The magnetization is plotted as a function of h/J± for J = 0.5Jj_ and for the values 
of J y /J± given in the legend. The (J/ J±, J x / J±) and (qi = J±/J,<%2 = Jx/J) phase diagrams 
at M z = 0.25 are displayed in (b) and (c), respectively. 



B. Magnetization and phase diagrams 

The magnetization M z calculated using BMFT is displayed as a function of magnetic 
field in Fig. |4ji for three sets of coupling values also used in Ref.— . Because all the states 
are gapped in zero field, the magnetization remains zero for all fields smaller than the 
zero-field gap. This is clearly observed in Fig. 3^,. Our curves for M z compare very 
well qualitatively with existing numerical results^. For certain coupling values a plateau 
appears at M z = 0.25 which means that on average half the spins are aligned along the 
magnetic field. Remember that M z = 0.5 is the saturation value. The frustration stabilizes 
this state, and as J x increases the size of the plateau also increases. The (J/J±, J x /J±) 
and (ai,a2) phase diagrams in Figs. HJd and [3J: show the region where the plateau at 
M z — 0.25 appears. Here, ot\ = Jj_/J and oli = J x /J- This plateau region is larger than, 
and encompasses, the numerically calculated plateau region^. We find that the plateau 
appears as J x increases only when J± > 0.67J. Also, when J x >> J± the plateau 
disappears. The sharp boundaries in the M z = 0.25 phase diagrams only appear at T = 0. 
At finite temperature they are replaced by crossovers. Note that the (J/ J±, J x /J±) phase 
diagram is symmetric with respect to the diagonal because the Hamiltonian is symmetric 
under exchanging the J and J x termsi. 

The effect of temperature on the magnetization, and on the magnetization plateau is 
illustrated in Fig. [5^. As temperature increases the size of the plateau decreases, and 
eventually disappears altogether at high enough temperature. Because of thermal excita- 
tions, spins originally locked in the gapped state at zero temperature become available for 
alignment along the magnetic field. This leads to the linear increase of magnetization in 
the low-field regime for T — 0.3J±/kB', Fig. [5^,. At T = there is a critical field at which 
the magnetization reaches the saturation value M s = 1/2, but for nonzero temperatures 
the critical behaviour is replaced by a crossover regime in a way similar to what was found 
for the Heisenberg chain and two- leg ladder in the absence of frustration 7 . 
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FIG. 5: An instantaneous configuration where on average the magnetization per site M z — 0.25. 
For simplicity all the spins on one chain are ferromagnetically arranged while on the other chain 
they are antiferromagnetically arranged. 

C. Spin bond parameters 

The bond-mean-field parameters are shown in Figs. [6p and [6^ for J = 0.5Ji and 
J x = 0.4 Jx. At the plateau, statistically half of the spins have been aligned with the 
magnetic field because M z — 0.25. Surprisingly, the bond-mean-field parameters also 
adopt half of their maximum value. The maximum value of Q and P' is their value when 
J_L = h = 0. The maximum value of P is 0.5. For simplicity imagine an instantaneous 
configuration where all the spins on one chain align along one direction, whereas the spin 
on the other chain which are antiferromagnetically arranged. This is illustrated in Fig. [5] 
In this case, the average magnetization per site is M z = 0.25, and exactly half of the 
rung bonds, one of the chains, and half of the diagonal bonds are antiferromagnetically 
arranged. The bond parameters represent the average AF correlations in the directions 
along which they are calculated, so each of them adopts half of its maximum value in zero 
field. The bond parameters vanish in the strong field regime, where the state labeling 
becomes obsolete and all three free energies equal. The bond parameters vanish at the 
same field where the magnetization reaches saturation signaling in this way a quantum 
(zero-T) phase transition^. 

When T > 0, thermal fluctuations spoil this transition. A crossover behaviour replaces 
the zero-T critical behaviour exactly in the same way as in the Heisenberg chain and 
non- frustrated two- leg ladder-. At high enough temperature the plateaus in the bond 
parameters disappear; Figs. [5p and[5p. 

D. Uniform spin susceptibility 

For convenience, in this section the unit of energy is J rather than Jj_. The magnetic 
susceptibility \ — ^jf- is shown as a function of field for three different couplings in 
Fig. [7K at temperature ksT/ J = 0.01J, and for three different temperatures in Fig. Wp 
for oti = J±/J = 1 and a-2 = J x /J = 0.5. For a% = 1 and a.2 — 0, \ is characterized 
by two peaks only, and is zero for low fields because of the energy gap. As a 2 increases 
X shows more structure and a gap. The gap in \ is due to the plateau at M z = 0.25. 
Thermal fluctuations round-off the peaks and close the gaps in \ m agreement with the 
disappearance of the plateau at high enough temperatures. 

E. Interpretation of the results 

The magnetization plateau and its behaviour in general can be explained using the 
chemical potential and energy band filling of the JW fermions in the same way as in Ref i 
The excitation energies £jyi and are plotted in figure [5^, as functions of wavenumber 
k for two different values of field h, and for J = 0.5J± and J x = QAJ±. When h = 0, the 
chemical potential of the JW fermions is zero; so only the lower-energy bands En% and 
Epj4 arc filled. A glance at Eq. (fT2"]) reveals that since two out of four bands are filled it 
is clear that M z — 0. As the field increases, the chemical potential of the JW fermions 
increases. For fields h smaller than the energy gap, M z remains zero. For the present set 
of exchange coupling constants, h must be greater than 0.57J^, which is the value of the 
energy gap in the absence of a magnetic field, in order to get a nonzero magnetization; 
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FIG. 6: The magnetization in (a) and the spin bond parameters in (b) and (c) are plotted as 
functions of magnetic field for three different temperatures; ksT — O.OlJx, fcgT = 0.1Jj_ and 
fcsT = 0.3J±. The exchange coupling constants are J = 0.5Jj_ and J x = 0AJ±. 




FIG. 7: The susceptibility \ ~ ls plotted as a function of h/J for (a) three coupling values 
at UbT I J = 0.01, and (b) three temperatures with ai = 1 and ai = 0.5. 



i.e., for fields greater than this threshold value the population of the higher-energy bands 
increases as field increases, leading to an increasing magnetization M z . At h — l.Uj_ the 
state of the system becomes the R-type state. The excitation energies E^i and Er 2 are 
plotted in Fig. [SJd as functions of k for two different values of field. The energy bands in 
the R-type state are separated by a gap for certain coupling values whereas in the other 
two states they are not. From h = l.Uj_ to h = 1.9J± three of the bands are populated 
and the fourth is empty so that there is a plateau at M z — 0.25. Further increasing 
the field increases the population of the fourth band until it becomes completely full at 
h = 2.42 Jj_. At this point M z = 0.5, and further increasing the field can no longer increase 
the population of the JW fermions. This results in the saturation of the magnetization 
M z . 
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FIG. 8: (a) The energies Emi and Em2 are plotted as functions of wavenumber k for h = 0J± 
and h — 0.5 J±. (b) The energies Em and Er2 are plotted as functions of k for h = 1.1 Jl and 
h = 2.45 J±. Here, J = 0.5 J±, and J x = 0.4 J±. 



IV. CONCLUSION 



In this work we analyzed the interplay between the quantum criticalities induced by 
frustration and magnetic field in the two-leg antiferromagnctic Hcisenberg ladder. When 
frustration is introduced in the presence of a magnetic field a plateau appears at half the 
saturation magnetization. The plateau in the magnetization appears because of the onset 
of an energy gap between two (excitation) energy bands. As frustration increases this 
plateau increases in size indicating that this phase is stabilized by the frustration. Our 
results agree qualitatively well with existing numerical data of Sakai and Okazaki. In the 
presence of frustration, the system is characterized by three spin bond parameters, one 
along the chains, one along the rungs and one along the diagonals. These parameters 
vanish in the strong field limit signaling the occurrence of a zero-T phase transition. 
The latter give place to a crossover behaviour when temperature becomes nonzero. The 
magnetization plateau shrinks as temperature increases, and disappears in the high-T 
regime. The spin-bond parameters are also characterized by a plateau that disappears in 
the high-T regime. 
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